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^ I Abstract 

Assuming that Quantum Einstein Gravity (QEG) is the correct theory of gravity 
Oh! on all length scales we use analytical results from nonperturbative renormalization 

. group (RG) equations as well as experimental input in order to characterize the 

special RG trajectory of QEG which is realized in Nature and to determine its 
^ ■ parameters. On this trajectory, we identify a regime of scales where gravitational 

I physics is well described by classical General Relativity. Strong renormalization 

effects occur at both larger and smaller momentum scales. The latter lead to a 
growth of Newton's constant at large distances. We argue that this effect becomes 
visible at the scale of galaxies and could provide a solution to the astrophysical 
missing mass problem which does not require any dark matter. We show that an 
extremely weak power law running of Newton's constant leads to flat galaxy rotation 
curves similar to those observed in Nature. Furthermore, a possible resolution of 
the cosmological constant problem is proposed by noting that all RG trajectories 
admitting a long classical regime automatically give rise to a small cosmological 
constant. 



1 Introduction 



During the past few years, in the hght of a series of investigations [1-9], it appeared 
increasingly likely that Quantum Einstein Gravity (QEG), the quantum field theory of 
gravity whose underlying degrees of freedom are those of the spacetime metric, can be 
defined nonpcrturbativcly as a fundamental, "asymptotically safe" [10] theory. By defini- 
tion, its bare action is given by a non-Gaussian renormalization group (RG) fixed point. 
In the framework of the effective average action [11-13] a suitable fixed point is known to 
exist in the Einstein-Hilbert truncation of theory space [1,3,6] and a higher-derivative 
generahzation [5] thereof. Detailed analyses of the reliability of this approximation [3-5] 
and a conceptually independent investigation [14] suggest that the fixed point should 
indeed exist in the exact theory, implying its nonperturbative renormalizability. 

The general picture regarding the RG behavior of QEG as it has emerged so far 
points towards a certain analogy between QEG and non-Abelian Yang-Mills theories. 
Quantum Chromo-Dynamics (QCD) say. For example, hke the Yang-Mills couphng 
constant, the running Newton constant G — G{k) is an asymptotically free coupling, it 
vanishes in the ultraviolet (UV), i.e. when the typical momentum scale k becomes large. 
In QCD the realm of asymptotic freedom, probed in deep inelastic scattering processes, 
for instance, is realized for momenta k larger than the mass scale Aqcd which is induced 
dynamically by dimensional transmutation. In QEG the analogous role is played by the 
Planck mass mpi. It delimits the asymptotic scaling region towards the infrared (IR). For 
k S> mpi the RG fiow is well described by its linearization about the non-Gaussian fixed 
point [4]. Both in QCD and QEG simple local truncations of the running Wilsonian action 
(efi^cctivc average action) are sufficient above Aqcd and mpi, respectively. However, as the 
scale k approaches Aqcd or mpi from above, many complicated, typically nonlocal terms 
arc generated in the effective action [8,15,16]. In fact, in the IR, strong renormalization 
effects arc to be expected because gauge (diffeomorphism) invariance leads to a massless 
excitation, the gluon (graviton), implying potential IR divergences which the RG ffow 
must cure in a dynamical way. Because of the enormous algebraic complexity of the 
corresponding fiow equations it is extremely difficult to explore the RG fiow of QCD or 
QEG in the IR, far below the UV scaling regime, by purely analytical methods. In QCD 
lattice techniques can be used to study the IR sector, but despite recent progress on 
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dynamical triangulations [17, 18] there exists no comparable tool for gravity yet. 

In QCD we have another source of information about its small momentum or large 
distance regime. If we take it for granted that QCD is the correct theory we can ex- 
ploit the available experimental data on the strong interaction, interpret them within this 
theory, and thus obtain information about the quantum dynamics of QCD, in particular 
its nonperturbative IR sector, from the purely phenomenological side. An example to 
which we shall come back in a moment are the non-relativistic quark-antiquark poten- 
tials extracted from quarkonium data (and confirmed on the lattice). They suggest that 
nonperturbative IR effects modify the classical Coulomb term by adding a confinement 
potential to it which increases (linearly) with distance: 

V(r)^-- + Kr. (1) 
r 

Here a and the string tension k are constants [19]. 

In this paper we are going to apply a similar "phenomenological" strategy to gravity. 
Under the assumption that QEG is the correct theory of gravity on all distance scales, 
we try to describe and characterize the distinguished RG trajectory which is realized in 
Nature as completely as possible. We use both observational input and the available 
analytical RG studies. 

We shall start from the flow equations in the Einstein-Hilbert truncation, determine 
which type of its RG trajectories the one realized in Nature belongs to, identify a regime 
on it where standard General Relativity is valid, and finally argue that this regime does 
not extend to arbitrarily large distances. In fact, for k smaller than the momenta typical 
of the regime of standard gravity, the truncation predicts a strong increase of G{k) with 
decreasing k, which, at a certain critical value of k, becomes infinite even. The diverging 
behavior is clearly an artifact of an insufficient truncation, but we shall see that the growth 
of Newton's constant with the distance can be understood on general grounds as due to 
a potential IR singularity. It is the main hypothesis of the present paper that a "tamed" 
form of this nonperturbative IR growth of G{k) is a genuine feature of exact QEG.^ 

The problem of the missing mass or "dark matter" is one of the most puzzling 
mysteries of modern astrophysics and cosmology [24]. It has been known for a long time 

^ Using different methods or models, IR quantum gravity effects have also been studied in refs. [20-23]. 
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that the luminous matter contained in a galaxy, for instance, does not provide enough 
mass to explain the gravitational pull the galaxy exerts on "test masses" in its vicinity. 
Typically their rotation curves v{r), the orbital velocity as a function of the distance, are 
almost flat at large distances rather than fall off according to Kepler's law [25]. Similar, 
but even stronger mass discrepancies are observed on all larger distance scales and in 
particular in cosmology. The recent high-redshift supernova and CMBR data show very 
impressively that the known forms of baryonic matter account only for a small percentage 
of the matter in the Universe. A possible way out is the assumption that the missing mass 
is due to some sort of "dark matter" which would manifest itself only by its gravitational 
effects. However, as to yet it has not been possible to convincingly identify any dark 
matter candidate, and so it might be worthwhile to think about alternatives. 

It is a very intriguing idea that the apparent mass discrepancy is not due to an 
unknown form of matter we have not discovered yet but rather indicates that we are using 
the wrong theory of gravity, Newton's law in the non-relativistic and General Relativity 
in the relativistic case. In fact, Milgrom [26] has developed a phenomenologically very 
successful non-relativistic theory, called MOdified Newtonian Dynamics or "MOND", 
which explains many properties of galaxies, in particular their rotation curves, in a unified 
way without invoking any dark matter. In its version where gravity (rather than inertia) 
is modified, a point mass M produces the potential 

0(r) = + ^JaoGM ln(r) (2) 

where G and Oq are constants. The second term on the RHS of Q is responsible for 
the fiat, non-Keplerian rotation curves at large distances. So far no wholly satisfactory 
relativistic extension of MOND is known. 

Also the relativistic theory proposed by Mannheim [27] where the Lagrangian is the 
square of the Weyl tensor tries to explain the rotation curves as due to a non-Newtonian 
force. The corresponding potential is of the form 

GM ^ 

(p[r) = \- nr. (3) 

r 

The resulting rotation curves do not become fiat but still seem to be in accord with the 
observations. 
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For a detailed discussion of other attempts at modifying gravity at astrophysical 
distances and a comprehensive hst of references we refer to [28]. A possible relation to 
quintessence has been speculated about in [29]. 

In the present paper we are going to explore the idea that the IR quantum effects 
of QEG, in particular the growth of G at large distances, induces a modified Newtonian 
potential similar to © or ©. If so, one can perhaps solve the missing mass problem in a 
very elegant and "minimal" manner by simply quantizing the fields which are known to 
exist anyhow, without having to introduce "dark matter" on an ad hoc basis. ^ 

It is particularly intriguing that the potentials we would like to derive within QEG 
are strikingly similar to the nonperturbative quark-antiquark potentials generated by 
(quenched) QCD. In particular eqs. ((T)) and © are mathematically identical, describing 
"linear confinement", while the MOND potential increases slightly more slowly at large 
distances. In view of the many similarities between QCD and QEG it is hard to believe 
that this should be a mere coincidence. 

The purpose of the present paper is to learn as much as possible about the gravi- 
tational RG trajectory Nature has chosen and to investigate the possibility that the IR 
renormalization effects of QEG are "at work" at galactic and cosmological scales. 

The remaining sections of the paper are organized as follows. In Section |21 we dis- 
cuss the Einstein-Hilbert truncation of theory space with an emphasis on the strong IR 
renormalization effects it gives rise to. In Sections El and |3J we analyze the RG trajec- 
tory realized in Nature, first using mostly analytical results (Section ^ and then also 
phenomenological input (Section Hj). In Section we employ a plausible model of the 
trajectory in the deep IR to demonstrate that an extremely tiny variation of Newton's 
constant would explain the observed fiat rotation curves. The results are summarized in 
Section [HI 



^See refs. [22,23] for a related analysis within a perturbatively renormalizable higher derivative gravity. 
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2 Towards the infrared with the 
Einstein— Hilbert truncation 



2.1 Structure of the RG flow 

In this subsection we discuss some properties of the Einstein-Hilbert truncation, in 
particular the classification of its RG trajectories. The emphasis will be on their behavior 
in the IR. We refer to [1] and [4] for further details. 

Our basic tool is the effective average action Fj. [gfj.u\ , a free energy functional which 
depends on the metric and a momentum scale k with the interpretation of a variable 
infrared cutoff. The action Ffe is similar to the ordinary effective action F which it ap- 
proaches for A; ^ 0. The main difference is that the path integral defining F^ extends only 
over quantum fiuctuations with covariant momenta > k^. The modes with < k"^ 
are given a momentum dependent (mass)^ oc Rk{p^) and are suppressed therefore. As a 
result, Ffc describes the dynamics of metrics averaged over spacetime volumes of linear 
dimension k~^. The functional F^ [g^u] gives rise to an effective field theory valid near 
the scale k. Hence, when evaluated at tree level, F^ correctly describes all quantum grav- 
itational phenomena, including all loop effects, provided the typical momentum scales 
involved are all of order k. 

Considered a function of k, F^ describes a RG trajectory in the space of all action 
functionals. The trajectory can be obtained by solving an exact functional RG equation. 
In practice one has to resort to approximations. Nonperturbative approximate solutions 
can be obtained by truncating the space of action functionals, i. e. by projecting the RG 
flow onto a finite-dimensional subspace which encapsulates the essential physics. 

The "Einstein-Hilbert truncation" , for instance, approximates F^ by a linear combi- 
nation of the monomials J^/g' R and J^/g^- Their prefactors contain the running Newton 
constant G{k) and the running cosmological constant A(/c). Their /c-dependence is gov- 
erned by a system of two coupled ordinary differential equations. 

The flow equations resulting from the Einstein-Hilbert truncation are most conve- 
niently written down in terms of the dimensionless "couplings" g{k) = k'^~'^ G{k) and 
\{k) = A{k)/k'^ where d is the dimensionality of spacetime. Parameterizing the RG tra- 
jectories by the "RG time" t = lak the coupled system of differential equations for g and 



5 



A reads dtX = (3^, dtg = f3g, where the /3-functions are given by 

X 2did+ 1) $i/2(-2A) - 8 d $i/2(0) - d{d + 1) r/N ^/^{-^X) 
Pg{X,g) = {d-2 + r]^) g 
Here r]-^, the anomalous dimension of the operator Jy/g' R, has the representation 

VN{g,X) = — — . 

l-gB2{X) 

The functions -Bi(A) and B2{X) are defined by 



(4) 



(5) 



d{d + 1) $^/2-i(-2A) -6d{d- 1) ^l^i-^X) 



d/2\ 



4d$i/2-i(0)-24<|.2/,(0) 



d/2\ 



(6) 



52(A) = -I (47r)i-'^/2 + 1) $i/2-i(-2A) - 6 - 1) $ 



-2A) 



The above expressions contain the "threshold functions" $^ and They are given by 



r(n) 



dz 



-1 



Z] — z 



(7) 



and a similar formula for $^ without the i?^°'''-term. In fact, R^^^ is a dimensionless version 
of the cutoff function i?^, i. e. R^ (p^) oc P (j9^/ A;^) . Eq. ((Tj) shows that ^^{w) becomes 
singular for w ^ — 1. (For all admissible cutoffs, z + R^^\z) assumes its minimum value 
1 at z = and increases monotonically for z > 0.) If A > 0, the $'s in the /3-functions 
are evaluated at negative arguments w = — 2A. As a result, the /3-functions diverge 
for A 1/2 and the RG equations define a flow on a half-plane only: — oo < 5^ < oo, 
-00 < A < 1/2. 

This point becomes particularly clear if one uses a sharp cutoff [4]. Then the <l>'s 
either display a pole at w = —1, 

1 1 1 



^ Y{n) p- 1 (1 + 



(8) 
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or, in the special case p = 1, they have a logarithmic singularity at w 

^l{w) = -T{n)-^ ln(l +w) + ipn- 



-1: 



(9) 



The constants ipn = ^n(O) parameterize the residual cutoff scheme dependence which 
is still present after having opted for a sharp cutoff. In numerical calculations we shall 
take them equal to the corresponding $;'j(0)-value of a smooth exponential cutoff^, but 
their precise value has no influence on the qualitative features of the RG flow [4]. The 
corresponding $'s are constant for the sharp cutoff: = 6pi/T{n + 1). 

From now on we continue the discussion in = 4 dimensions. Then, with the sharp 
cutoff, the coupled RG equations assume the following form: 



dtX = - (2 - //n) A - 



TT 



dtg 



(2 + r/N) g 
6tt + 5 g 



5 ln(l - 2 A) - (^2 + I r?N 



+ 5 ln(l-2A) -^1 + 6 



(10a) 
(10b) 
(10c) 



1 - 2A 

Solving the system (jlOj) numerically [4] we obtain the phase portrait shown in Fig. 
n] The RG flow is dominated by two flxed points ((?*, A,,,): a Gaussian flxed point (GFP) 
at (7* = A* = 0, and a non-Gaussian fixed point (NGFP) with > and A* > 0. 
There are three classes of trajectories emanating from the NGFP: trajectories of Type 
la and Ilia run towards negative and positive cosmological constants, respectively, and 
the single trajectory of Type Ila ( "separatrix" ) hits the GFP for A; ^ 0. The short- 
distance properties of QEG are governed by the NGFP; for — > oo, in Fig. ^ all RG 
trajectories on the half-plane g > run into this point. The conjectured nonperturbative 
renormalizability of QEG is due to the NGFP: if it is present in the full RG equations, it 
can be used to construct a microscopic quantum theory of gravity by taking the limit of 
infinite UV cutoff along one of the trajectories running into the NGFP, thus being sure that 
the theory does not develop uncontrolled singularities at high energies [10]. By definition, 
QEG is the theory whose bare action 5* equals the fixed point action limfc^oo T/c [gfiu] ■ 

The trajectories of Type Ilia have an important property which is not resolved in 
Fig. ^ Within the Einstein-Hilbert approximation they cannot be continued all the way 

■^For this purpose we employ the exponential cutoff with "shape parameter" s = 1. In = 4, the only 
(fi's we need are (pi = C(2) and Lp2 = 2C(3). See ref. [4] for a detailed discussion of the sharp cutoff. 
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Figure 1: RG flow on the (?-A-plane. The arrows point in the direc- 
tion of decreasing values of k. (From ref. [4].) 



down to the infrared {k = 0) but rather terminate at a finite scale fcterm > 0. At this scale 
they hit the singular boundary A = 1/2 where the /3-functions diverge. As a result, the 
flow equations cannot be integrated beyond this point. The value of /cterm depends on the 
trajectory considered. 

In ref. [4] the behavior of g and A close to the boundary was studied in detail. 
The aspect which is most interesting for the present discussion is the following. As the 
trajectory gets close to the boundary, A approaches 1/2 from below. In this domain the 
anomalous dimension ()10c|) is dominated by its pole term: 

on + 5g 1 — 2A 

Obviously r/N \ —oo for A 1/2, and eventually r^N = —oo at the boundary. This 
behavior has a dramatic consequence for the (dimensionful) Newton constant. Since 
dtG = r/N G, the large and negative anomalous dimension causes G to grow very strongly 
when k approaches fcterm from above. This behavior is sketched schematically in Fig. |21 

At moderately large scales k, well below the NGFP regime, G is approximately 
constant. As k is lowered towards kterm, G{k) starts growing because of the pole in oc 
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1/ (1 — 2 A), and finally, at k = /cterm, it develops a vertical tangent, (dG/dk) (/cterm) = 
— oo. The cosmological constant is finite at the termination point: A(A;term) = ^torm/2- 





G(k) 




Glab — 






Fig 


kterm kj^i, 

are 2: Schematic behavior of G{k) for 
trajectories of Type Ilia. 



By fine-tuning the param- 
eters of the trajectory the scale 
fcterm Can be made as small as 
we like. 

Since it happens only ve- 
ry close to A = 1/2, the diver- 
gence at kterm is not visible on 
the scale of Fig. [TJ (Note also 
that g and G are related by a 
decreasing factor of k'^.) 



The phenomenon of trajec- 
tories which terminate at a finite scale is not special to gravity, it occurs also in truncated 
fiow equations of theories which are understood much better. Typically it is a symptom 
which indicates that the truncation used becomes insufficient at small k. In QCD, for 
instance, thanks to asymptotic freedom, simple local truncations are sufficient in the UV, 
but a reliable description in the IR requires many complicated (nonlocal) terms in the 
truncation ansatz. Thus the conclusion is that for trajectories of Type Ilia the Einstein- 
Hilbert truncation is reliable only well above ktevm- It is to be expected, though, that in 
an improved truncation those trajectories can be continued to = 0. 

We believe that while the Type Ilia trajectories of the Einstein-Hilbert truncation 
become unreliable very close to /ctcrnu their prediction of a growing G{k) for decreasing k 
in the IR is actually correct. The function G{k) obtained from the differential equations 
(fTUI) should be reliable, at least at a qualitative level, as long as A <^ 1. For special 
trajectories the IR growth of G{k) sets in at extremely small scales k only. Later on we 
shall argue on the basis of a gravitational "RG improvement" [30-35] that this IR growth 
is responsible for the non-Keplerian rotation curves observed in spiral galaxies. 

The other trajectories with g > 0, the Types la and Ila, do not terminate at a finite 
scale. The analysis of ref. [4] suggests that they are reliably described by the Einstein- 
Hilbert truncation all the way down to k = 0. 
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2.2 What drives the IR renormaUzation? 



Next we discuss a simple physical argument which sheds light on the dynamical 
origin of the expected strong IR effects. As we shall see, they are due to an "instability 
driven renormalization" , a phenomenon well known from many other physical systems 
[16,36], spontaneous symmetry breaking being the prime example. 

For an arbitrary set of fields 0, and in a slightly symbolic notation^, the exact RG 
equation for the effective average action reads [11,12] 



9trJ0] = iTr 



2 



(r[')[0] + i?,) 'dtRi\ . (12) 



It contains the fully dressed effective propagator (r^.^-* + fZ^) ^ where r^^"* denotes the 
Hessian of and Rk the cutoff operator. It is instructive to rewrite ()12|) in the form 

5*r.[</>] = i-^lndet(rf [</>]+ i?,) (13) 

where the derivative D / Dt acts on the /c-dependence of Rk only. The RHS of ()13|) repre- 
sents a "/3-functional" which summarizes all the infinitely many ordinary /3-functions in 
a compact way. Obviously its essential ingredient is a kind of a "one-loop determinant" . 
It differs from that of a standard one-loop calculation by the presence of the cutoff term 
Rk and by the use of the dressed inverse propagator F^.^-* rather than the classical S^'^\ 
More importantly, in the standard situation one expands the classical action S[(f)\ about 
its minimum 0min(a;), in which case the one-loop effective action lndet(S'*^^^[0min]) sums 
up the zero-point energies of the small stable oscillations about 0min(a;). On the RHS of 
flT^ . instead, (pi^x) is a prescribed external field, the argument of F^ on the LHS. It can 
be changed freely, eq. (fT^ holds for all 4>{x), so that (p{x) is not in general a stationary 
point of Ffc. 

Thus we may conclude that the basic physical mechanism which drives the RG flow 
is that of quantum fluctuations on arbitrary off-shell backgrounds 0. They determine the 
/3-functional, and depending on how "violent" those fluctuations are, the RG running is 
weaker or stronger. 

It is helpful to consider two extreme cases. Let us first assume that, for a certain fixed 



^We ignore possible complications due to gauge invariance. They are inessential for the present 
discussion. 
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(f){x), the operator is positive definite; then +Rk is positive, too,^ and the quadratic 
action governing the fluctuations 6(j){x) about (j){x) given by Jd^x 6(f) (Ff ) + e) 6(f) IS a 
positive quadratic form. In this case the computation of the /3-functionaI amounts to 
summing up the zero-point energies of small stable fluctuations which would not grow 
unboundedly. The "fluctuation induced" renormalizations of the parameters in F^. which 
they give rise to are comparatively weak. 

(2) 

If, on the other hand, F^, [0] has one or several negative eigenvalues fi < with 
< oo, then F^.^-* + Rf. = F^^-* + k'^ is positive only in presence of the IR regulator, for 
k'^ > Without the IR regulator there is a real physical instability. Within the linear 
approximation the fluctuation modes grow unboundedly; beyond the linear approximation 
they would try to "condense" in order to turn (f> into a stable ground state. Following a 
well behaved RG trajectory one stays in the regime k'^ > where F^^"* + k^ is positive. 
However, when k"^ approaches |yu| from above, the lowest eigenvalue of F^^"* + fc^ gets 
very close to zero, and the RG flow is strongly affected by the presence of the nearby 
singularity. The effective propagator (F^^-* + Rk) ^ becomes very large, and typically this 
leads to an enormous growth of the (standard) /3-functions when k^ \ They give rise 
to comparatively strong "instability induced" renormalizations. We shall see in a moment 
that the IR effects of QEG are precisely of this type. 

In order to flnd the RG flow on the full theory space the above stability analysis 
and the "summation of zero-point energies" has to be performed for inflnitely many 
different backgrounds (j){x)] they are needed in order to "project out" all the possible fleld 
monomials which constitute the functional F^. On a truncated theory space just a few 
0's might be sufficient. 

In order to illustrate the relationship between the (ordinary) /3-functions and the 
instability presented by F^.^-* let us look at a scalar model (on flat spacetime) in a simple 
truncation: 

Tk[4>] = jd^x [ld,^d>'4>+lm\k)(f)' + ^^X{k)(f)^y (14) 

Here (f> denotes a real, 22^symmetric scalar fleld, and the truncation ansatz ()14|) retains 
only a running mass and (/)^-coupling. In a momentum basis where —d^d^ = > we 

^In order to capture the essence of the argument it is sufficient to use a mass-type cutoff [12] for 
which Rk = k'^ and R° {p^/k"^) = 1 for all p^. 
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have 

r^') ^p^ + m\k) + \{k)(P^. (15) 
Always assuming that A > 0, we see that T^^^ is positive if > 0; but when < it 



k 

k2 



can become negative for small enough. Of course, the negative eigenvalue for = 0, 
for example, indicates that the fluctuations want to grow, to "condense", and thus to 
shift the field from the "false vacuum" to the true one. By the mechanism discussed 
above, this gives rise to strong instability induced renormalizations. In fact, the standard 
/3-functions for and A can be found by inserting (jl5j) into (jl2j) . taking two and four 
derivatives with respect to 0, respectively, and then setting = in order to project out 
dtm"^ and dtX. As a result, the /3-functions are given by p-integrals over (powers of) the 
propagator 

p"^ + m^{k) + k"^ \ (16) 

In the symmetric phase (m^ > 0) this (euclidean!) propagator has no pole, and the 
resulting /3-functions are relatively small. In the broken phase (m^ < 0), however, there 
is a pole at = —m{kY — k"^ provided k"^ is small enough: k"^ < \m{kY\. For k"^ \ |m(fc)^| 
the /3-functions become large and there are strong instability induced renormalizations. 

In a reliable truncation, a physically realistic RG trajectory in the spontaneously 
broken regime will not hit the singularity at fc^ = |m(fc)^|, but rather make m(/c) run in 
precisely such a way that |m(A;)^| is always smaller than k"^. This requires that 

-m(A;)2 oc k'^. (17) 

This strong instability induced mass renormalization is necessary in order to evolve an 
originally W^-shaped symmetry breaking classical potential into an effective potential 
which is convex and has a fiat bottom. (See [12] for a detailed discussion of this point.) 

Unfortunately the two-parameter truncation ()14|) is too rudimentary for a reliable 
description of the broken phase. Its RG trajectories actually do run into the singularity. 
They terminate at a finite scale /cterm with k^^^^ = \m{kterm)^\ at which the /3-functions 
diverge. Instead, if one allows for an arbitrary running potential f/fc(0), containing in- 
finitely many couplings, all trajectories can be continued to = 0, and for A; \ one 
finds indeed the quadratic mass renormalization p7|) [12]. 
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Let us return to gravity now where corresponds to the metric. In the Einstein- 
Hilbert truncation it suffices to insert the metric corresponding to a sphere S^(r) of arbi- 
trary radius r into the flow equation^ in order to disentangle the contributions from the 
two invariants / d'^Xy/g' oc and J d'^Xy/g' R oc r^. Thus we may think of the Einstein- 
Hilbert flow as being a manifestation of the dynamics of graviton fluctuations on S^(r). 
This family of backgrounds, labeled by r, is "off-shell" in the sense that r is completely 
arbitrary and not fixed by Einstein's equation in terms of A. 

It is convenient to decompose the fluctuation hf^i, on the sphere into irreducible (TT, 
TL, ■ ■ ■ ) components [3] and to expand the irreducible pieces in terms of the corresponding 
spherical harmonics. For /i^,^ in the transverse-traceless (TT) sector, say, the operator 
Tfc + -Rfc equals, up to a positive constant, 

-D^ + 8r-^ + -2A{k) (18) 

with D"^ = g^'^ D^Dy the covariant Laplacian acting on TT tensors. The spectrum of 
— D^, denoted {p^}, is discrete and positive. Obviously (fTH|l is a positive operator if the 
cosmological constant is negative. In this case there are only stable, bounded oscillations, 
leading to a mild fluctuation induced renormalization. This is precisely what we observe in 
the IR of the Type la trajectories: there is virtually no non-canonical parameter running 
below k = rripi. The situation is very different for A > where, for k"^ sufficiently small, 
flTH|) has negative eigenvalues, i. e. unstable eigenmodes. In fact, expanding the RHS of the 
flow equation to orders and the resulting /3-functions are given by traces (spectral 
sums) containing the propagator [1] 

[p^ + e - 2 A{k)]~\ (19) 

The crucial point is that the propagator ()19|) can have a pole when A (A;) is too large and 
positive. It occurs for A(A;) > k'^/2, or equivalently A(A;) > 1/2, at p"^ = 2 A{k) — k"^. Upon 
performing the p^-sum this pole is seen to be responsible for the terms oc 1/ (1 — 2 A) and 
ln(l — 2 A) in the /3-functions which become singular at A = 1/2. The allowed part of 
the g-X-plane (A < 1/2) shown in Fig. ^ corresponds to the situation k"^ > 2A{k) where 
the singularity is avoided thanks to the large regulator mass. When k"^ approaches 2A{k) 

^We stress, however, that the /3-functions do not depend on the choice of background used for 
projecting onto the various invariants [1,3,12]. 
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from above the /3-functions become large and strong renormalizations set in, driven by 
the modes which would go unstable'' at /c^ = 2 A. 

In this respect the situation is completely analogous to the scalar theory discussed 
above: Its symmetric phase (m^ > 0) corresponds to gravity with A < 0; in this case 
all fluctuation modes are stable and only small renormalization effects occur. Conversely, 
in the broken phase (m^ < 0) and in gravity with A > 0, there are modes which are 
unstable in absence of the IR regulator. They lead to strong IR renormalization effects 
for k"^ \ \m{ky \ and k"^ \ 2A{k), respectively. The gravitational Type la (Type Ilia) 
trajectories are analogous to those of the symmetric (broken) phase of the scalar model. 

As for the behavior of the RG trajectories near the boundary the crucial question 
is whether, when k is lowered, A{k) decreases at least as fast as k"^ or more slowly. In the 
first case the trajectory would never reach the singularity, while it does so in the second. 
In the Einstein-Hilbert approximation the trajectories of Type Ilia indeed belong to the 
second case; since A{k) does not decrease fast enough the RG trajectory runs into the 
pole at a certain fcterm where fcterm — 2A(A;term)- 

The termination of certain trajectories is not specific to gravity. We saw that it 
happens also in the scalar model if we use the over-simplified truncation ()14|) . This simple 
ansatz has similar limitations as the Einstein-Hilbert truncation. In the scalar case the 
cure to the problem of terminating trajectories is known [12]: If one uses a more general 
truncation, allowing for a non-polynomial Uk{(f)), the RG trajectories never reach the 
singularity and extend to A; = 0, with strong renormalizations, however, in particular the 
quadratic running (fT7|) . In a certain sense the Einstein-Hilbert truncation has a similar 
status as a polynomial truncation for Uk{4>)'- it is not general enough to be reliable down 
to A; = for a positive cosmological constant or in the broken phase, respectively. While 
there are computationally manageable truncations of sufficient generality in the scalar 
case it is not known which truncations would allow for a reliable continuation of the Type 
Ilia trajectories below fcterm- They are likely to contain nonlocal invariants [8,16,39] which 
are hard to handle analytically. 



^From the propagator H19|l it is obvious that the smaUer the eigenvalue of a fluctuation mode, 
the higher is the scale k'^ at which this particular mode starts contributing significantly, and the more 
important is its impact on the RG flow. A galaxy-size fluctuation is more important than a solar system- 
size fluctuation, for example. 
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In view of the scalar analogy it is a plausible and very intriguing speculation that, 
for — i> 0, an improved gravitational truncation has a similar impact on the RG flow as 
it has in the scalar case. There the most important renormalization effect is the running 
of the mass: — m(A;)^ oc k"^. If gravity avoids the singularity in an analogous fashion the 
cosmological constant would run proportional to fc^, 

A{k) = \'^e (20) 

with a constant A^^ < 1/2. In dimensionless units fl20j) reads X{k) = Aj,^, i.e. A^^ is a 
fixed point of the A-evolution. If the behavior (j2(Jj) is actually realized, the renormalized 
cosmological constant observed at very large distances, A{k — ^ 0), vanishes regardless of 
its bare value. Clearly this would have an important impact on the cosmological constant 
problem [40]. 

3 The RG trajectory realized in Nature: 
the Einstein— Hilbert domain 

How can we find out which one of the RG trajectories shown in Fig. ^ is realized 
in Nature? As in every quantum field theory, one has to experimentally determine the 
value of appropriate "renormalized" quantities.^ In Quantum Electrodynamics (QED), 
for instance, one measures the electron's charge and mass in a large-distance experiment, 
thus fixing e{k) and me{k) at = 0. In QCD the point A; = is inaccessible, both 
theoretically and experimentally, so one uses a "renormalization point" at a higher scale 
k > 0. In QEG the situation is similar. In the extreme infrared {k 0) we have 
neither theoretically reliable predictions nor precise experimental determinations of the 
gravitational couplings g and A. 



^Above we considered pure gravity, while the experimental data include renormalization effects due to 
matter fields. Since in this paper we are interested in order of magnitude estimates only, and we anyhow 
do not know the exact matter field content of Nature, we assume that the inclusion of matter does not 
change the general qualitative features of pure gravity. As for the nonperturbative rcnormalizability it is 
known that there exist matter systems with this property and that they are "generic" in a sense [7]. 



15 



3.1 Exploiting experimental information 

We know that all gravitational phenomena at distance scales ranging from terres- 
trial experiments to solar system measurements are well described by standard General 
Relativity (GR). Since this theory is based upon the Einstein-Hilbert action with con- 
stant values of G and A, we can conclude that the RG evolution of those parameters for 
k between the related typical mass scales (1 meter) ~^ and (1 astronomical unit)~^, say, 
is negligibly small. In this "GR regime" the renormalization group flow is essentially 
the canonical one, i.e. g{k) oc k"^ and X{k) oc l/k"^ which follows from g{k) = k'^G{k), 
X{k) = A{k)/k'^ when G,A = const. The corresponding RG trajectories are the hyper- 
bolas g oc 1/X depicted in Fig. |21 During the /c-interval defining the GR regime the true 
RG trajectory realized in Nature must be very close to one of those hyperbolas. 

Which class does this trajectory 
belong to? Recent CMBR and high 
redshift supernova data show that the 
present Universe is in a state of accel- 
erated expansion which, in a Fried- 
mann-Robertson- Walker framework, 
can be explained by an nonzero pos- 
itive cosmological constant. In the 
RG context this should mean that 
A(A; ^ Hoj > since the relevant 
scale is set by Hq, the present Hub- 
ble parameter. Among the trajecto- 
ries of the Einstein-Hilbert trunca- 
tion only those of Type Ilia and Illb 
run towards positive A's for \ 0. 
Trajectories of Type Illb correspond to a negative G and are excluded therefore. Thus 
the RG trajectory realized in Nature, as long as it remains in the domain of validity of 
the Einstein-Hilbert truncation, belongs to Type Ilia. 

We saw that trajectories of Type Ilia cannot be continued below a certain /cterm; 
and in the present paper we are going to argue that fcterm is roughly of the order of typical 
galaxy scales. As a result, the Einstein-Hilbert truncation is probably insufficient to de- 
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Figure 3: The canonical RG flow cor- 
responding to constant values of the 
dimensionful parameters G and A. 
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Figure 4: The Type Ilia trajectory realized in Nature and the separatrix. The 



scribe the (continuation of the) Type Ilia trajectory realized in Nature at the cosmological 
scale k ^ Ho where A was actually measured. Nevertheless it seems to be clear that, for k 
large enough, the true RG trajectory selected by Nature is a Einstein-Hilbert trajectory 
of Type Ilia. The reason is that the other alternatives. Type la and Type Ila, can be 
computed reliably down to /c = 0, and they do not give rise to a positive cosmological 
constant in the infrared. 

For our picture to be correct, the prospective Type Ilia trajectory must contain a 
sufficiently long "GR regime" where it runs on top of one of the hyperbolas of Fig. |21 The 
situation is sketched qualitatively in Fig. EJ The Type Ilia trajectory spirals out of the 
NGFP, approaches the separatrix, runs almost parallel to it for a while, then "turns left" 
near the GFP, and finally runs towards the singularity at A = 1/2. After the turning point 
where c^^A = 0, but before it gets too close to A = 1/2, this trajectory is an almost perfect 
hyperbola of the canonical RG fiow. In Fig. E] we indicate the latter by the dashed line 
which, between the points Pi and P2, is indistinguishable from the Type Ilia trajectory. 
It is this segment between Pi and P2 which can be identified with the realm of classical 
GR. As we shall see in a moment, the variation of G and A along the true Type Ilia 
trajectory is unmeasurably small between Pi and P2. 

Which one of the infinitely many Type Ilia trajectories did Nature pick? We can 
answer this question if, from experiments or astrophysical observations, we know a single 



dashed line is a trajectory of the canonical RG fiow. 
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point {g, A) the trajectory passes through. Let us assume we measure the (dimensionful) 
Newton constant and cosmological constant in a "laboratory" with a typical linear di- 
mension of the order l/fciab- We interpret the result of the measurements as the running 
couplings evaluated at this scale: G(fciab), ^{huh)- Knowing those two values, as well as 
the pertinent "laboratory" scale fciab, we can compute the dimensionless couplings: 

^(Mab) = kl^G{h^^), A(A;iab) = A(A;iab)/A;L- (21) 

The pair ((7(/ciab), A(A;iab)) uniquely fixes a trajectory in the Einstein-Hilbert approxima- 
tion. If one uses a more general truncation, further parameters need to be measured, of 
course. The first one of the eqs. (j2H) can be rewritten in the following suggestive form: 

9{kub) = (A;iab/mpi)' = {ipi/KlT- (22) 

Here we defined the Planck length and mass in the usual way in terms of the measured 
Newton constant G[k\a,h) according to ip\ = mpj^ = \J G{k\ab) ■ 

Newton's constant has been measured at length scales k^J^ ranging from the size 
of terrestrial experiments to solar system dimensions. Within the errors the result has 
always been the same: G'(fciab) = G\ab = 6.67- 10^^^ m^ kg^^ s^^. We can now calculate 
^(^lab) according to (jUJ). At the typical scale of a terrestrial laboratory one finds 

^7(fciab) ~ 10-'° forfc-^ = lm (23) 

while at the solar system scale of 1 astronomical unit , 

g{h^^,)^lQ-^^ ioik^J^ = lA\J. (24) 

In any case g{kiah) is an extremely small number for any fciab in the GR regime, g{k\a.h) ^ 
1. Its precise value will not matter in the following; for the sake of clarity we shall use the 
example k^J^ = 1 meter and g{kia.\^) ^ 10^™ for numerical illustration. Throughout the 
discussion the length scale k^l is assumed to lie in the GR regime, ranging from terrestrial 
to solar system distances. 

The determination of the associated A(A;iab) is difficult; in fact, rather than at "lab- 
oratory" scales 1 m ■ ■ ■ 1 AU, A was actually measured at cosmological distance scales. For 
a first qualitative discussion the following estimate is sufficient, however. According to the 
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effective vacuum Einstein equation at the scale fciab, a cosmological constant of magnitude 
■^(^lab) leads to the spacetime whose radius of curvature is of the order 

We know that in absence of matter, at k^J^ = 1 m, say, spacetime (when observed with 
a "microscope" of resolution k^J^) is flat with a very high precision, i. e. that r^. is much 
larger than the size of the "laboratory": Vc ^ k^J^. As a consequence, A(/i;iab) must be 
very small compared to unity: 

A(fciab) < 1. (25) 

This means in particular that, at fciab and in the entire GR regime, the trajectory realized 
in Nature is still very far away from the dangerous singularity at A = 1/2 where the 
Einstein-Hilbert approximation breaks down. 



3.2 Approximate RG flow near the GFP 

Since (?, A ^ 1 in the GR regime we may neglect higher order terms g^,gX, • ■ ■ in 
the eqs. (|lUj) and obtain the following flow equation linearized about the GFP: 

dtX = -2 X + 1^29/7^ (26a) 
dt9 = 2g. (26b) 

In the "linear regime" where ()26p is valid, only the cosmological constant shows a non- 
canonical running, while the dimensionful Newton constant does not evolve in the ap- 
proximation ()26b|) . 

In the next to leading order G runs according to dtG = rj^G with t]-^ = —bg 
proportional to g. In any cutoff scheme 6 is a positive constant of order unity, b = 
(24 — v^i) /Svr for the sharp cutoff. Hence 

-77N ^ 10-™- ■■ 10-^1 (27) 

RG regime 

The smallness of these numbers explains the success of standard General Relativity based 
upon the approximation ?7n = 0, and it confirms our interpretation of the segment between 
the points Pi and P2 in Fig. |3] as the realm of classical gravity. 
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To the right of the point P2 in Fig. ^ at scales k lower than those of the GR regime, 
the growth of G{k) due to the infrared instability sets in. 

Let us return to the linearized flow equations (|26|). They are applicable whenever the 
trajectory is close to the GFP {g, A ^ 1), not only in the GR regime. We may use them to 



which it switches from decreasing to increasing values of A. (See Fig. El) Setting dfX = 
in (|26a|) we obtain 



The constant is cutoff scheme-, i. e. i?*^°''-dependent, but it is of order unity for any 
cutoff. Therefore 



Eqs. (jlHl), are vahd provided qt, \t ^ 1- Later on we shall see that this is actually 
the case in Nature. 

After the trajectory has passed the turning point, g keeps decreasing and A increases. 
As a result, the second term on the RHS of ()26ap . {p2g/TT, gradually becomes negligible 
compared to the first one, —2 A, so that the flow equation becomes the canonical one. 
This marks the beginning of the GR regime at Pi. 

It is easy to solve the coupled differential equations (j26j) exactly. They allow for two 
free constants of integration which we fix by requiring gikr) = Qt and \ {kT) = \t- By 
definition, fcy is the scale at which the trajectory passes through the turning point. The 
solution reads 



derive a relation between the coordinates (^gx, At) of the trajectory's turning point T at 



At = (v^2/27r) gr- 



(28) 



Xr/gr = 0{1). 



(29) 




(30a) 




(30b) 



The corresponding running of the dimensionful parameters is given by 



G(k) = -pr = const 



(31a) 




(31b) 
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Since the hnear regime contains the GR regime, we may identify the constant in 
with Giab = nipf. This entails the important relation k'^ = gxrapi, or 

kr = 



9t mpi. 

We observe that a small qt will lead to a large hierarchy mpi ^ k^- 
We can use fl32|) in order to eliminate kx from ()30|): 

2 



(32) 



mpi 



mpiv 



1 + 



2 4 



For later use we note that for any k in the linear regime 

1 



G{k) A{k) = g{k) X{k) = - grX 



T 



k 

kx 



(33a) 
(33b) 

(34) 



Looking at eqs. ()31b .b) we see that, while G does not run at all in the linear regime, 
the scale dependence of A is entirely due to the factor [l + (/c/fc^-)^] . Once k has become 
much smaller than fc-r, after the trajectory has "turned left" , this factor approaches unity, 
and A effectively stops to run. By definition, this happens at Pi, the starting point of the 
GR regime. 

Let us make this statement more precise. Denoting by ki the scale at which the 
trajectory passes through Pi, the requirement is that {ki/kr)'^ <^ 1. We quantify the 
precision with which the fc^-term is negligible in the GR regime by means of an exponent 
i>. In terms of i>, we define ki by 

ki/kr = 10-". (35) 

As a result, {k/kx)^ is smaller than 10"^'^ for all scales in the GR regime {k < ki). A 
value such as v = 1 should be sufficient in practice. It makes sure that in the GR regime 
A is constant with a precision better than 0.01 %. 

For k any scale in the GR regime we obtain from ()33p. in very good approximation. 



g{k) = [k/mpiY 

X{k) = IqtXt [rripi/ky. 



(36a) 
(36b) 
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Similarly pip yields the following constant values of the dimensionful quantities: 

G{k) = Giab (37a) 

A{k) = \\Tkl = \K{kT). (37b) 

Remarkably, in the GR regime, A differs from its value at the turning point precisely by 
a factor of 1/2. As a trivial consequence, 

(GA) = i (GA) (38a) 

GR regime turning point 

and (jH^ reduces to 

G{k)k{k)=g{k)\{k) = \gTXT (38b) 

We shall come back to this relationship later on. 

Next we derive various estimates for g and A at kx and ki. If we evaluate (jHUjl at 
k = ki, neglecting the fc'^-term in ()30b|) . and use (jHSjl we find^ 

g{k,) =gT 10"'"' (39a) 

A(fci) = At10+2^ (39b) 

Obviously Pi has a g- (A-) coordinate which is smaller (larger) than the corresponding 
coordinate of T by a factor 10~'^'^ (10"'"^'^), as it should be according to Fig. 0] Now 
we exploit eq. (j^^ . Neglecting factors of order 1 we have gx = Ay which yields when 
combined with 

g{k,)=X{k,)lO-^''. (40) 

Since A increases along the trajectory we know that A^ < X{ki) < A(A;iab) where 
k\ah is any "laboratory" scale in the GR regime. Therefore, as a consequence of the 
experimental result (j2Sl), 

At < A(fci) < A(A;iab) < 1- (41) 
Since Xt and gx are almost equal, gx is small, too: 

5(T ~ At < 1. (42) 



'in order~of-magnitude equations such as (|39|) we suppress inessential factors of order unity. 
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According to ()4U|) . g{ki) is even smaller than A(fci). Hence, with X{ki) ^ 1 from ()4ip . it 
follows that 

(?(A;i)<l, A(A;i)<l. (43) 

Eqs. (|42p and (j43p show that for the RG trajectory realized in Nature the points T 
and Pi are located at an extremely short distance to the GFP. The trajectory starts at 
the NGFP with coordinates g^.,X^. = (9(0.1) [3,5]. Then it follows the separatrix until, 
at very tiny values of g and A, it gets ultimately driven away from the GFP along its 
unstable A-direction. In pictorial terms we can say that the trajectory is squeezed deeply 
into the wedge formed by the separatrix and the g = 0-axis. As a consequence, it spends 
a very long RG time near the GFP because the /3-functions are small there. In this sense 
the RG trajectory which Nature has selected is highly non-generic or "unnatural". It 
requires a precise fine-tuning of the initial conditions, to be posed infinitesimally close to 
the NGFP. 

3.3 Existence of a GR regime and the 
cosmological constant problem 

Why did Nature pick a trajectory which gets so "unnaturally" close to the GFP? 
Why not, for example, one of those plotted in Fig. ^which always keep a distance of order 
unity to the GFP and require no special fine-tuning? Of course questions of this kind 
cannot be answered within QEG, for the same reason one cannot compute the electron's 
charge or mass in QED. 

However, it is fairly easy to show that a Universe based upon one of the generic 
trajectories would look very different from the one we know. The main difference is that, 
along a generic trajectory, no sufficiently long GR regime would exist where classical 
gravity makes sense at all. According to our previous discussion, the GR regime is located 
in between a regime with strong UV renormalization effects (spiraling around the NGFP 
related to asymptotic safety) and, most probably, a second regime with a significant 
running of the parameters in the IR. For classical GR to be applicable the UV and IR 
regimes must be well separated. For a generic trajectory this is not the case, however. 
The generic Einstein-Hilbert trajectories computed in [8] leave the UV regime at /c ~ mpi 
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and soon after they terminate at a fcterm not much smaller than mpi; there is no GR regime 
which would last for a few orders of magnitude at least. 

The basic mechanism which allows for the emergence of a GR regime is to fine-tune 
the trajectory in such a way that it spends a lot of RG time near the GFP. In this manner 
the onset of the IR regime, within the Einstein-Hilbert truncation characterized by the 
value of kterm, gets enormously delayed, fcterm being much smaller than mpi. If classical 
GR is correct up to a length scale L, the trajectory must be such that its /c^rm larger 
than L since at k^J.^ the IR effects are likely to become visible. 

A quantitative estimate of fcterm can be obtained as follows. In the linear regime 
the RG flow is explicitly given by eqs. ()30p. Once k is sufficiently low we leave this 
regime and the full nonlinear Einstein-Hilbert flow equations have to be used. At even 
smaller scales, the truncation breaks down and we should switch to a more general one. 
Outside the linear regime, the /3-functions are no longer small, hence the trajectory has a 
comparatively high "speed" there. As a result, it takes the trajectory much longer to go 
from T to the boundary of the linear, i. e. GR regime than from there to a point close to 
A = 1/2. Therefore we may use eq. ()36b|) for \{k) in the GR regime in order to derive an 
estimate for /cterm- We identify fctcrm with the scale where, according to ()36b|) . the value 
A = 1/2 is reached: Qt^t ('^Pi/^term)^ = 1- In this rough approximation, 

kterm = \/ Qt^t ?^pi (44a) 
= {ip2/2TiY''^ grmpi (44b) 



or, in terms of length scales. 



^term " / ■, (^5) 



This equation shows explicitly that by making qt ~ At small, k^^^ can be made as large 
as we hke. 

Thus we have demonstrated that only a "non-generic" trajectory with an "unnatu- 
rally" small qt ^ Xt ^ does give rise to a long GR regime comprising many orders of 
magnitude. 

In we saw that the turning point is passed at kx = y/gj^mpy. Ignoring the 
C(l)-factor (y92/2vr) in ()44b|) this entails fcterm = JWrkx- As a result, there exists an 
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exactly symmetric "double hierarchy" among the three mass scales fcterm, ^t, and mp\: 

qt < 1, = < 1- (46) 



krp m-pi 

Therefore, on a logarithmic scale, kf is precisely in the middle between mpi and /cterm- 
Thanks to the smallness of qt it is many orders of magnitude away from either end. 

The emergence of a long regime where gravity is essentially classical is one of the 
benefits we get from the unnaturalness of the trajectory chosen by Nature. Another one is 
that, in this classical regime, the cosmological constant is automatically small. Inserting 
(I32|) into ()37b|) we obtain the following A in the GR regime: 

A (A;) = \gT^Ti^'p\ = const. (47) 

Again thanks to the smallness of qt and Ay, the cosmological constant is much smaller 
than m|j. Up to a factor of order unity. 



A 



= gl « 1- (4^ 

GR regime 



Thus we may conclude that the very fine-tuning which gives rise to a long GR regime at 
the same time implies a large hierarchy between the A in this GR regime and m|j, which 
often is considered its "natural" value. 

This observation provides a solution to the "cosmological constant problem" by 
realizing that it is actually part of a much more general naturalness problem. Rather 
than "Why is A so small?" the new question is "Why does gravity behave classically over 
such a long interval of scales?" . 

The basic reason for this connection is very simple. Denoting the cosmological 
constant in the GR regime by A(/ciab) = Aiab we have \{k) = Ai^h/k'^ there. This X{k) 
approaches unity so that the IR renormalization effects become strong once k is of the 



order vAiab • Hence, roughly. 



"^term 



(49) 



which shows that fcterm is small if, and only if, Aiab is small. 

This is precisely what one finds by numerically solving the flow equations: There 
do not exist any Type Ilia trajectories which, on the one hand, admit a long classical 
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regime, and on the other hand, have a large cosmological constant. As a consequence, 
the smallness of the cosmological constant poses no naturalness problem beyond the one 
related to the very existence of a classical regime in the Universe. 

Strictly speaking this resolution of the cosmological constant problem is only a par- 
tial one for the following reasons. (1) We analyzed the flow equations of pure gravity 
only, but we believe that the inclusion of matter fields (forming symmetry breaking con- 
densates, etc.) will not change the general picture. Since anyhow only massless particles 
contribute to the /3-functions for A; ^ 0, it is hardly possible that the matter fields de- 
stroy the IR renormalizations near A = 1/2. But if they survive the estimate ()49|) remains 
intact, again implying that a long GR regime requires a small Aiab- (2) Our argument 
refers to Aiab rather than the cosmologically relevant K{Hq). Because of the IR effects, 
Aiab and A{Ho) differ in principle, but we shall see later on that in Nature this difference 
is small compared to the notorious 120 orders of magnitude one has to cope with. 

In fact, the upper boundary of the linear regime is roughly the Planck scale since, 
according to (j^ . g{k) and X{k) are of the order of g^ and A* slightly below k = mpi. 
Between this scale and the turning point the change of A is enormous. From ()31b|) with 
(1221) we obtain A(mpi) = jA{kT)/gT ^ A(A;t)- In the next section we shall see 

that gr ^ 10~^°, whence A(A;t) ~ 10^^^° A(mpi). 

3.4 Is the Hubble scale within the classical regime? 

The observational (CMBR, supernova, etc.) data, when interpreted within standard 
cosmology, show that the present vacuum energy density of the Universe is very close to 
the critical one, implying that A is of the order of Hq. (The general relationship is 
A = 3QaHq, and the data favor I^a ~ 0.7). 

We can now ask whether, if A is as large as Hq, the Hubble scale is still within the 
GR regime. If we interpret Hq as a "laboratory" value, Aiab ~ Hq, the estimate (jl^ 
yields 

^term ^ Hq. (50) 

This is a very remarkable and intriguing result: On the RG trajectory Nature has picked, 
the Hubble scale k = Hq is precisely at the boundary of the GR regime. At distances small 
compared to the Hubble length £h = Hq'^ classical GR is a good approximation, but on 
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length scales C ^ the IR renormalization effects become important so that its use is 
questionable there. 

Since k = Hq is just at the boundary of the GR regime it is plausible to assume 
that the cosmological values G{Hq) and A(ifo) do not differ too much from Giab and Aiab, 
respectively. (Later on we shall see that the difference should not be more than one or 
two orders of magnitude.) If so, we are in a position to completely fix the parameters 
of the RG trajectory because we can fit the measured A-value to the trajectories of the 
linearized flow. Inserting k = Hq, A(/c) = Hq, and G{k) = mpf into ()38b|) we find 

lgTXT = 9{k)X{k)= {Ho/mpi)\ (51) 

Since qt and At are approximately equal, and since the "experimental" value for Hq/tupi 
is known to be about 10"^'', eq. ()49|) allows for an explicit determination of qt and Xt'- 

qt^Xt^ Ho/mpi ^ 10"^°. (52) 

This number is indeed extremely small but, consistently with the picture drawn so far, 
still larger than 5'(/i;iab) ~ 10^™. 

With the "double hierarchy" comprises 30 orders of magnitude between 
fcterm and kT, and between kr and mpi, respectively. Combining eq. ()44a|) for A;term with 
we rediscover that /cterm ~ Hq, and using eq. we obtain the scale at the turning 
point: 

kr ^ ^/Hq mpi ^ 10-^°mpi. (53) 

The associated length scale is 

k^^ ^ 10^° ipi ^ 10^^ cm. (54) 

It is very exciting that this is a truly macroscopic length scale, very far away from the 
Planck regime which is usually thought to be the realm of the quantum gravitational UV 
renormalization effects. Note that the "turning left" of the trajectory, despite the large 
value of (j54p . is still an effect of the "UV type" in the sense that it has nothing to do 
with the IR phenomena occurring when A ~ 1/2. It happens so extremely late because 
of the enormous fine-tuning of the trajectory. 



27 



Allowing a margin of one order of magnitude (i^ = 1), say, ()54|) implies a beginning 
of the GR regime at 




10^2 cm. 



(55) 



At length scales ^ = k larger than this value A is constant with high precision. But 
for £ ^ A;f^ it has a strong scale dependence oc k'^ given by eq. ()31b|) . Since \Tk\j2 = 
Aiab ~ Hq we may rewrite this equation in the form 



Should we then expect to see strong violations of GR at the millimeter or micrometer 
scale? The answer is probably no. While the variation of A (A;) is indeed very strong, 
between k~^ = 1 /im and k~^ = 1 nm, say, it changes by 12 orders of magnitude, the 
absolute value of A is proportional to Hq which is extremely tiny according to all standards 
of terrestrial experiments, of course. Because of the practical problems with measuring 
A at non-cosmological distances the running of A could possibly remain undetected over 
many orders of magnitude. 

By combining with we obtain g and A at the Hubble scale: 



This approximation neglects the IR effects close to Hq, but should be correct as far as the 
orders of magnitude are concerned. 

4 The RG trajectory realized in Nature: the deep IR 

4.1 Leaving the GR regime 

Denoting the scales corresponding to the boundary points of the GR regime. Pi and 
P2, by ki and ^2, respectively, the general picture of the trajectory which Nature has 
chosen can be summarized as follows. For = 00 it starts infinitesimally close to the 
NGFP, runs very close to the separatrix until, at k^, it turns left in the very last moment 
before hitting the GFP, at Qti^t ^ 1, then enters the GR regime at the scale ki and 
leaves it at ^2- In the entire GR regime G and A are constant, and \{k) ^ 1. 



A{k)^Hl l+{k/kTY 



(56) 



^7(//o) -10 



120 



(57) 
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At scales k < k2 immediately below the GR regime there exists a regime where the 
renormalization effects become appreciable, but are still weak enough for the Einstein- 
Hilbert truncation to be reliable. The mechanism discussed in Section |21 causes G{k) 
to increase with decreasing k. The dimensionless cosmological constant A is still much 
smaller than 1/2 there. 

At even smaller scales, A approaches 1/2, and the Einstein-Hilbert truncation 
becomes unreliable as it would yield an unbounded growth of G{k) and r/N —oo, 
dG/dk — oo at a nonzero fcterm- 

It is our main hypothesis that exact QEG or a sufficiently general truncation "tames" 
this divergent behavior and leads to a bounded growth of G{k), i.e. a finite ?7n, along a 
RG trajectory which passes smoothly through fcterm and can be continued to A; = 0. We 
believe that those IR renormalization effects lead to observable physical effects, and that 
the scale at which they occur is approximately given by the value of kterm obtained from 
the Einstein-Hilbert truncation. 

The truncations needed for a reliable ab initio computation of the IR effects within 
QEG are much more complicated than those used so far. They contain many more 
couplings beyond g and A. The resulting RG trajectory can be visualized as a curve in 
a high dimensional "theory space" containing the ^f-A-plane as a subspace. In general 
one has to distinguish the projection of the exact trajectory onto the g-X-plane and the 
trajectory implied by the Einstein-Hilbert truncation whose theory space consists of this 
plane only. As long as this truncation is reliable the two curves in the (/-A-plane do not 
differ much, but the deviations become large when the trajectory runs out of the domain 
of validity of the truncation. For the RG trajectory realized in Nature this is the case 
near the question mark in Fig. |3J close to A = 1/2. 

Clearly it would be desirable to employ such improved truncations which allow 
for a continuation of the Type Ilia trajectories to = 0, but because of the enormous 
mathematical complexity of the resulting flow equations this seems to be out of reach 
using the presently available computational techniques. The situation is similar to QCD 
where the analogous nonperturbative IR phenomena such as color confinement or chiral 
symmetry breaking are notoriously difficult to deal with analytically. 

It is of crucial importance to find out at which scale the IR effects become strong. 
In principle /cterm can be determined by measuring G{kiah) and A(A;iab) in the "laboratory" 
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and integrating the Einstein-Hilbert flow equations towards smaller k until one reaches 
A = 1/2 at fcterm- However, in practice this does not work because it is very hard to 
measure A in a non-cosmological "laboratory" . 

A first hint about the value of fcterm comes from the argument in Subsection 13.41 
Given the measured value of A, fcterm is of the order of the Hubble scale. We take this 
as a strong indication suggesting that at least in cosmology the IR effects are "at work" 
already. However, in order to assess their relevance at galactic scales, a more precise 
estimate is needed. 

Since we have no analytical tools yet to investigate those IR effects we can resort 
to the following phenomenological strategy for obtaining information about the RG tra- 
jectory: We make an ansatz for the trajectory, work out its consequences, and compare 
them to the observational data.^° The so-called IR fixed point model developed in [32,33] 
is an attempt in precisely this direction. In the next subsection we describe its status in 
the context of the present discussion. 

4.2 The IR fixed point model 

The IR fixed point model [32, 33] is a cosmology based upon the hypothesis that, 
for k 0, the projected trajectory {g{k), \{k)) runs into a non-Gaussian IR fixed point 
[g^, X]^) , different from the UV fixed point it emanates from. This assumption implies 
that in the deep IR the dimensionful parameters run as 

G{k)=g'J'/e, A{k)=Xfe. (58) 

The behavior of A (A;) is the same as in (j2(Jj) which was motivated by the analogy with the 
scalar theory, and G{k) oc l/Zc^ is a plausible ansatz for the growth of Newton's constant. 
Since /3g = (2 + t^n) g holds true in general, but not necessarily with the r^N of (0), the 
anomalous dimension at the fixed point is ?7n = —2. 

Note that while the IR renormalization according to (j58p has a strong effect on G{k) 
and A{k) separately, it leaves the product G{k) A{k) invariant: A(A;) decreases at the same 

^'^In this context it is important to understand that the RG trajectory is a universal object which 
is not specific to any particular system: rfe[(7^,y] defines an effective field theory for all systems whose 
typical scale is of the order of k. The specific features of the system under consideration enter only when 
it comes to identifying k in terms of dynamical or geometrical data, see [31,37,38]. 
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rate G{k) increases towards the IR. 

In a homogeneous and isotropic Universe the relevant cutoff scale is = ^/t to 
leading order^'^ [31,32]. Here t is the cosmological time and ^ is a constant of order 
unity. This "cutoff identification" turns G{k) and A(A;) into functions of time: G{t) = 
G{k = ^/t), and similarly for A. Replacing the constants G and A in the classical 
Einstein equation with G{t) and A{t) one obtains a "RG improved" field equation whose 
Robert son- Walker-type solutions can be studied [31,32]. For the fc-dependence one 
finds an accelerating Universe with a scale factor a{t) oc t^/^. The fixed point not only 
forces the Universe to enter an epoch of accelerated expansion for t ^ oo, it also explains 
without any fine-tuning of parameters why in the late Universe the matter energy density 
equals approximately the vacuum energy density, thus providing a natural solution to 
the "cosmic coincidence puzzle". Confronting the infrared fixed point model with the 
observational data (supernovae, compact radio sources, CMBR, etc.) it performs as well 
as the best-fit Friedmann model. (See ref. [33] for further details.) 

The fixed point solution to the improved Einstein equation exists only if X]^ and C, 
are such that Aj,^^ ^ = 8/3, implying that Aj,^ has to be of order unity since ^ = 0{1). 
Furthermore, is found to be of the order [G'(to)/Giab] {Ho/mpif. Here G(to) denotes 
the present Newton constant at cosmological scales. According to the analysis of ref. [33], 
the observational data imply that the ratio G'(to)/G'iab can comprise at most one or two 
orders of magnitude (in any case a number by far smaller than 120, say). Therefore, using 
Ho/mpi ^ 10~^°, we find that the model is consistent only provided 

gf'^lO-''', Xf^O{l). (59) 

This model is based upon the assumption that the IR effects lead to the formation 
of a fixed point into which the trajectory is attracted for k —>■ 0. The solubility of the 
improved Einstein equation and the observational data then require that the projection 
of the fixed point onto the gf-A-plane is located at the coordinates ()59|) . The comparison 
with ()57|) shows that the point [g^^, \^^) is precisely in the region where the Einstein- 
Hilbert approximation breaks down and something new must happen. (Near the question 
mark in Fig. 0J) 

This consistency is an important success of the fixed point model. It describes the 
^^We consider only power law solutions for which f/t is proportional to the Hubble parameter H{t). 
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simplest possible behavior of the (projected) RG trajectory in the deep IR: the strong 
quantum corrections near A = 1/2 simply bring the running of g{k) and A(A;) to a complete 
standstill, at g^J^ and Aj,^, respectively. 

We can think of the cosmological time evolution as an evolution of the Universe along 
the RG trajectory with k corresponding to 1/t or H{t), which is essentially the same here. 
As long as the Universe stays in the GR regime the classical Friedmann equations are a 
valid description, but as soon as A gets close to 1/2 the fixed point cosmology takes over, 
the Universe starts accelerating, and the vacuum and matter energy densities become 
approximately equal {Qa ~ ^m)- These predictions are in remarkable agreement with 
what we learned about the Universe from the astrophysical data which became available 
during the past few years. (See refs. [32,33] for further details.) 

Even if the finer details of the IR fixed point model are perhaps not completely 
correct quantitatively it demonstrates that an increase (decrease) of G (A) near the Hub- 
ble scale would explain many of the observed features of the late Universe in an ex- 
tremely natural way. In a sense, the renormalization effects would mimic the presence of 
"quintessence" which was invented in order to explain these features. 

4.3 Where does the IR running set in? 

Motivated by the estimate fcterm ~ Hq from Subsection 13 . 41 and the phenomenological 
success of the IR fixed point model we believe that at the present Hubble scale the IR 
effects have an observable and qualitatively important impact on gravitational physics 
already. Now the crucial question is at which scale ki > Hq precisely the GR regime ends 
and the parameters start running. 

If we knew the exact RG trajectory we could set up the corresponding RG improved 
Friedmann equations, interpret the observational data measured in the late Universe 
within this framework, and deduce ki from the redshift of the epoch in which deviations 
from classical cosmology become visible. In ref. [33] this analysis was performed within 
the fixed point model but the statistical quality of the presently available data (on high 
redshift type la supernovae) is still too poor to allow for a precise determination of ki. It 
became fairly clear, however, that between k = ki and k = Hq Newton's constant cannot 
have changed by more than 1 or 2 orders of magnitude. 

Thus we are led to look for possible manifestations of the IR running in gravitation- 
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ally bound structures which are small on cosmological scales. Clearly the natural place 
to search for such phenomena are galactic systems. One of their most striking features, 
as to yet completely unexplained, is the discrepancy between their directly observable 
mass (due to "luminous matter" ) and the mass inferred from the observed motions. This 
mass discrepancy occurs in all types of galactic systems, from dwarf galaxies to super- 
clusters. Its magnitude ranges from a factor of a few at the kiloparsec scale to a factor 
of a few hundred at the megaparsec scale [24]. The most popular attempt at explaining 
this discrepancy is the dark matter hypothesis. 

In the present paper we propose that the mass discrepancy is actually not due to 
the presence of dark matter but rather to the IR renormalization effects. Taking this 
hypothesis seriously, we shall try to learn something about the RG trajectory by taking 
advantage of what is known about galaxies. There is a large amount of observational 
data [25] showing that the orbits of "test particles" moving in typical spiral galaxies 
cannot be explained by Newtonian gravity if the gravitational attraction is due to the 
luminous matter only. Usually the rotation curve v{r), the velocity on a circular orbit 
as a function of the distance to the center, becomes approximately constant for large 
distances r, in contradiction to Kepler's law. 

In the following sections we shall demonstrate, by means of an appropriate renor- 
malization group improvement, that such violations of Kepler's law are exactly what one 
would expect to occur if G{k) shows a certain scale dependence at galactic scales already. 

A first argument demonstrating that this scenario is not completely unreasonable is 
the following. A mass discrepancy of a factor of a few hundred at the cluster scale indicates 
that, very roughly, G at those scales should not differ from Giab by more than a similar 
factor of a few hundred. This is nicely consistent with what we got for G(to)/G'iab in the 
fixed point model [33] at somewhat larger scales. Given the many orders of magnitude 
we are dealing with here this consistency is certainly nontrivial. 

A comparatively small ratio Cgaiaxy/Giab oi 0(10) or perhaps O(IOO), and a similar 
ratio in cosmology, justifies our earlier use of the linearized flow equation down to k — Hq 
which neglects the running between the end of the GR regime and the Hubble scale. 

While a renormalization of G and A by a factor of 10 is extremely httle compared 
to the many orders of magnitude their values changed in the UV, it leads to signiflcant 
modifications of classical Newtonian gravity "in the sky" , at astrophysical scales. 
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4.4 Running G at galactic scales 

Henceforth we shall assume that the IR running starts somewhere between solar 
system and galactic scales. Concentrating on the running of Newton's constant^^ and, 
for technical simplicity, spherically symmetric model galaxies, the first step towards their 
RG improved dynamics is as follows. Given the function G = G{k) pertaining to the RG 
trajectory realized in Nature we convert its /^-dependence to a distance dependence by 
means of an appropriate cutoff identification [30]. As long as curvature effects are small, 
the obvious choice for spherically symmetric systems is k = ^/r. with ^ a constant of 
order unity [38]. We define the position dependent Newton constant as 

G{r) = G{k = i/r). (60) 

By the very construction of the effective average action [12], G{r) is the parameter which 
appears in the effective field equations for quantities averaged over a volume of linear 
dimension ^ r. 

Using G{r) we can try to RG improve the classical Newton potential —GiahM/r by 
substituting G\ah G{r). In the case of the UV renormalization effects [1] it is known 
that this procedure reproduces the results of the explicit perturbative calculation [41,42]. 
However, as we shall discuss in the next section the improved action approach predicts a 
"nonperturbative" large distance correction which is even more important than the one 
taken into account by replacing Giab G{r) in the potential. 

The most important question is which /c-dependence we should expect at a typical 
galactic scale of k'^ = 1 kpc ■ ■ ■ 100 kpc, say. As it will turn out, for the trajectory realized 
in Nature this scale is already outside the domain of validity of the nonlinear Einstein- 
Hilbert flow equations. Therefore, the best we can do is to make an ansatz for G{k) 
which reproduces the "phenomenology" of galaxies as well as possible, and at the same 
time provides (part of) a natural interpolation between the GR regime and the IR fixed 
point behavior. 

^•^As for the QCD analogy, a feature of QEG not shared by QCD is the long intermediate classical 
regime, as a result of which deviations from the 1/r-potential occur near fcterm ^ mpi only, while in 
QCD this happens at the scale Aqcd which is analogous to mpi in other respects. 

^^For small isolated systems such as an individual galaxy the cosmological constant plays no important 
role probably. (Cosmological effects will be neglected in our discussion of galaxies.) 
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Classical GR is characterized by a vanishing anomalous dimension, r^N = 0. Switch- 
ing on the renormalization effects, the theoretically simplest option is to assume a nonzero, 
but constant ?7n. Since dtG = rj^G this entails a power law G{k) oc k^^ . As G{k) is an 
increasing function of fc, t^n is negative. We shall set r/N = —Q therefore with a positive 
constant q. Thus we tentatively assume that 

G{k) oc l/A;" (61) 

for the range of fc-values which are relevant to the structure of galaxies. Note that 
the assumption of an approximately constant ?7n is not too restrictive since this range 
comprises only 2 or 3 orders of magnitude, from 1 kpc and 100 or 1000 kpc, say. The 
power law (j6ip leads to the radial dependence 

G{r) oc r^. (62) 

The consequences of this r-dependence we are going to explore in a moment. It will lead 
to phenomenologically acceptable rotation curves if g ^ 10^^. 

There exists an intriguing scenario which would give rise to an approximately con- 
stant ?7n during some /c-interval in a natural way. Let us assume that, at galactic scales, 
the RG trajectory approaches a further NGFP, different from the ones discussed above, 
spends some time in its vicinity, and is finally driven away from it. This fixed point 
A*, ■ ■ ■) must have the (^-coordinate = Q because then (3^ = (2 + r]j<f) g vanishes 
there without constraining the value of r]-^. (Every NGFP with (7* 7^ has r]^ = —2.) 
Hence, if 77=,, =77^(5'* = 0, A,,,, ■ ■ ■ ) = —q is nonzero, the smallness of the /3-functions near 
the NGFP implies that r]^{g{k), \{k), ■ ■ ■) is nonzero and approximately /c-independent 
there. 

A fixed point of this kind could only arise within truncations which are much more 
general than the ones used so far, as the result of a strongly nonperturbative dynamics. 
(Of course, the perturbative quantization of the Einstein-Hilbert action yields t^n = at 
g = 0, to all orders.) 

The projection of the new fixed point onto the ^f-A-plane of Fig. ^ is located exactly 
on the horizontal {g = 0) axis. What makes this scenario particularly plausible is that we 
know already that the RG trajectory realized in Nature does indeed get extremely close 
to the g = it would be necessary to make ?7n approximately /c-independent. 
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However, even if this additional NGFP should not exist as a strict fixed point it is still 
quite plausible that a behavior similar to (jHU) can prevail for a short RG time and this is 
all we need. 



5 Galaxy rotation curves 

In this section we apply the improved-action approach to the problem of the fiat 
galaxy rotation curves. For further details about the general approach we refer to [37], 
and to [38] for a detailed analysis of spherically symmetric, static "model galaxies" . 

The idea is to start from the classical Einstein-Hilbert action ^eh = /d^a^ \J~9 -^eh 
with the Lagrangian = (-^ ~ 2 A) / (IGvr G) and to promote G and A to scalar fields. 
This leads to the modified Einstein-Hilbert (mEH) action 

S.E„b.G,Al^^/d*.y3^{^-2A||}. (63) 

The resulting theory has certain features in common with Brans-Dicke theory; the main 
difference is that G{x) (and A(x)) is a prescribed "background field" rather than a Klein- 
Gordon scalar as usually. Upon adding a matter contribution the action (|63|) implies the 
modified Einstein equation^^ 

G^, = -A(x) g^, + 8n Gix) (T^, + AT^,) . (64) 

Here AT^j, is an additional contribution to the energy-momentum tensor due to the x- 
dependence of G: 

^T,u ^ ^ {D,D^ - g,. D') (65) 

The field equation ()64j) is mathematically consistent provided A{x) and G{x) satisfy a 
"consistency condition" which insures that the RHS of ()64|1 has a vanishing covariant 
divergence. 

In ref. [38] we analyzed the weak field, slow-motion approximation of this theory for 
a time-independent Newton constant G = G'(x) and A = 0. In this (modified) Newtonian 



^^In [37] and [38] a further contribution, was added to the energy-momentum tensor in order to 
describe the 4-momentum of the field G{x). Its form is not completely fixed by general principles. But 
as it does not affect the Newtonian limit [38] we set 9ni, = here. 
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limit the equation of motion for massive test particles has the usual form, x(t) = — V0, 
but the potential obeys a modified Poisson equation, 



V20 = 47rGpeff 



(66a) 



with the effective density 



.15 



PcS = P+{8nG) VW. 



(66b) 



In deriving it was assumed that T^j, describes pressureless dust of density p and that 
G(x) does not differ much from the constant G. We use the parameterization 



and assume that A/'(x) ^ 1. More precisely, the assumptions leading to the modified 
Newtonian limit are that the potential 0, the function A/", and typical (squared) velocities 
are much smaller than unity; all terms linear in these quantities are retained, but higher 
powers (0^, ■ ■ ■ ) and products of them {4>Af, ■ ■ ■ ) are neglected. (In the application to 
galaxies this is an excellent approximation.) Apart from the rest energy density p of 
the ordinary ( "baryonic" ) matter, the effective energy density Pefr contains the "vacuum" 
contribution 



which is entirely due to the position dependence of Newton's constant. Since it acts as a 
source for (p on exactly the same footing as p it mimics the presence of "dark matter". 

As the density (jUHj) itself contains a Laplacian V^, all solutions of the Newtonian 
field equation (j66|l have a very simple structure: 



Here (p is the solution to the standard Poisson equation = Av: G p containing only the 
ordinary matter density p. The simplicity and generality of this result is quite striking. 

Up to this point the discussion applies to an arbitrary prescribed position depen- 
dence of Newton's constant, not necessarily related to a RG trajectory. At least in the 

^^In this paper we neglect in pcff an inessential term M p relative to p. 




(67) 



(SttG) ' VW(x) = (SttG') Vg(x) 



(68) 



0(x) = 0(x) + |Ar(x). 



(69) 
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case of spherically symmetric systems the identification of the relevant geometric cutoff 
is fairly straightforward, oc 1/r, so that we may consider the function G{k) as the pri- 
mary input, implying G{r) = G{k = ^/r). Writing again G = G [1 + M] we assume that 
G{k) is such that A/" ^ 1. Then, to leading order, the potential for a point mass reads, 
according to 

^^r) = -^ + \U{r). (70) 

Several comments are in order here. 

(a) The reader might have expected to find a term —GMN'{r)/r on the RHS of (f7n|) 
resulting from Newton's potential 0n = —GM/r by the "improvement" G G{r). 
However, this term ^nA/" is of second order with respect to the small quantities we are 
expanding in. In the envisaged application to galaxies, for example, ^nA/" is completely 
negligible compared to the |A/'-term in (jTU)). 

(b) According to ()70|) . the renormalization effects generate a nonclassical force (per unit 
test mass) given by — A/''(r)/2 which adds to the classical 1/r^-term. This force is at- 
tractive if G{r) is an increasing function of r and G{k) a decreasing function of k. This 
is in accord with the intuitive picture of the antiscreening character of quantum grav- 
ity [1]: "Bare" masses get "dressed" by virtual gravitons whose gravitating energy and 
momentum cannot be shielded and lead to an additional gravitational pull on test masses 
therefore. 

(c) The solution (j7Up is not an approximation artifact. In [38] we constructed exact 
solutions of the full nonlinear modified Einstein equations (with JV not necessarily small) 
which imply (f7n|) in their respective Newtonian regime. Those exact solutions can be 
interpreted as a "deformation" of the Schwarzschild metric (M 7^ 0) or the Minkowski 
metric (M = 0) caused by the position dependence of G. The solutions related to the 
Minkowski metric are particularly noteworthy. They contain no ordinary matter (no point 
mass), but describe a curved spacetime, a kind of gravitational "soliton" which owes its 
existence entirely to the x-dependence of G. At the level of eq. (f7n|) they correspond to the 
M = 0-potential = | A/" which solves the modified Poisson equation if the contribution 
oc V^A/" is the only source term. In the picture where dark matter is replaced with a 
running of G this solution corresponds to a pure dark matter halo containing no baryonic 
matter (yet). The fully relativistic M = 0-solutions might be important in the early 
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stages of structure formation [38]. 

Let us make a simple model of a spherically symmetric "galaxy" . For an arbitrary 
density profile p = p{r) the solution of eq. reads 

0(r)= ^jdr'^^^ + \M{r) (71) 

r 

where Ai{r) = Air jdr' r' p{r') is the mass of the ordinary matter contained in a ball of 



radius r. We are interested in periodic, circular orbits of test particles in the potential 
fl71|) . Their velocity is given by v'^{r) = r(f)'{r) so that we obtain the rotation curve 

n, , GM(r) 1 d , 

v^{r) = — + 7:r — Af{r). 72 

r 2 ar 

We identify p with the density of the ordinary luminous matter and model the 
luminous core of the galaxy by a ball of radius Tq. The mass of the ordinary matter 
contained in the core is A^(ro) = Aio, the "bare" total mass of the galaxy. Since, by 
assumption, p = and hence A4{r) = Aio for r > ro, the potential outside the core is 
{/)(r) = —GAio/r + A/'(r)/2. We refer to the region r > tq as the "halo" of the model 
galaxy. 

As an example, let us adopt the power law G{k) oc which we motivated in Section 
m We assume that this fc-dependence starts inside the core of the galaxy (at r < vq) so 
that G{r) oc r'^ everywhere in the halo. For the modified Newtonian limit to be realized, 
the position dependence of G must be weak. Therefore we shall tentatively assume that 
the exponent q is very small (0 < g ^ 1); applying the model to real galaxies this will 
turn out to be the case actually. Thus, expanding to first order in g, r'^ = l + g ln(r) + ■ ■ ■ , 
we obtain G{r) = G[l+ J\f{r)] with 

7V(r) = q In(Kr) (73) 

where k is a constant. In principle the point G about which we linearize is arbitrary, but 
in the present context G = Giab is the natural choice. In the halo, eq. (fT3|l leads to a 
logarithmic modification of Newton's potential 

cPir) = -^^^ + UniKr). (74) 



39 



The corresponding rotation curve is 

2/ ^ GMq q 

V r = + -. 75 

r 2 

Remarkably, at large distances r oo the velocity approaches a constant Voo = \fqj2 . 
Obviously the rotation curve implied by the A;~'^-trajectory does indeed become flat at 
large distances - very much like those we observe in Nature. 

Typical measured values of range from 100 to 300 km/sec so that, in units of the 
speed of light, foo ~ 10~^. Thus, ignoring factors of order unity for a first estimate, we 
find that the data require an exponent of the order 

q ^ 10"^ (76) 

The smallness of this number justifies the linearization with respect to A/". It also implies 
that the variation of G inside a galaxy is extremely small. The relative variation of 
Newton's constant from some ri to r2 > ri is AG/G = q ln(r2/ri). As the radial extension 
of a halo comprises only 2 or 3 orders of magnitude the variation between the inner and 
the outer boundary of the halo is of the order AG/G ~ g, i. e. Newton's constant changes 
by one part in a million only. 

Including the core region, the complete rotation curve reads 

2, , GM(r) q 

= ^ + |. (77) 

The r-dependence of this velocity is in qualitative agreement with the observations. For 
realistic density profiles, Ai{r)/r is an increasing function for r < ro, and it decays as 
Aio/r for r > tq. As a result, f^(r) rises steeply at small r, then levels off, goes through 
a maximum at the boundary of the core, and finally approaches the plateau from above. 
Some galaxies indeed show a maximum after the steep rise, but typically it is not very 
pronounced, or is not visible at all. The prediction of (|75|) for the characteristic r-scale 
where the plateau starts is 2G M-^/q] at this radius the classical term G M-q/t and the 
nonclassical one, g/2, are exactly equal. With q = 10~^ and A4o = 10^^ Mq one obtains 
9 kpc, which is just the right order of magnitude. 

The above f^(r) is identical to the one obtained from standard Newtonian gravity 
by postulating dark matter with a density pdm oc 1/r^. We see that if G{k) oc k~'^ with 
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q ^ 10~^ no dark matter is needed. The resulting position dependence of G leads to an 
effective density pefr = P + <?/ (STrGr^) where the 1/r^-term, well known to be the source 
of a logarithmic potential, is present as an automatic consequence of the RG improved 
gravitational dynamics. 

We consider these results a very encouraging and promising indication pointing in 
the direction that quantum gravitational renormalization effects could be the origin of 
the plateaus in the observed galaxy rotation curves. If so, the underlying RG trajectory 
of QEG is characterized by an almost constant anomalous dimension rj-^ = —q ^ —10^^ 
for k in the range of galactic scales. 

Is the Einstein-Hilbert truncation sufficient to search for this trajectory? Unfortu- 
nately the answer is no. According to eq. (jlUcp . ?7n is proportional to g which is extremely 
tiny in the regime of interest, smaller than its solar system value 10"^^. In order to achieve 
a It^nI as large as 10~^, the smallness of g must be compensated by large IR enhancement 
factors. As a result, A should be extremely close to 1/2, in which case the RHS of pOc|) 
is dominated by the pole term: rj-^ ^ — [6g/7c) (1 — 2 A)~^. Assuming g ^ 10"^^ as a 
rough estimate, a g-value of 10~^ would require 1 — 2A ~ 10^®^. It is clear that when 
1 — 2 A is so small the Einstein-Hilbert trajectory is by far too close to its termination 
point to be a reliable approximation of the true one. Moreover, ?7n ((?(/;;), A(A;)) is not 
approximately fc-independent in this regime. Thus we must conclude that an improved 
truncation will be needed for an investigation of the conjectured RG behavior at galactic 
scales, in particular to search for the g^ = 0-fixed point speculated about at the end of 
Section HI 

It is clear that the above model of a galaxy is still quite simplistic and does not 
yet reproduce all phenomenological aspects of the mass, size, and angular momentum 
dependence of the rotation curves for different galactic systems. In particular Voo is a 
universal constant here and does not obey the empirical Tully-Fisher relation. Similarly, 
according to MOND, the need for a dark matter substitute actually does not arise at a 
critical distance but at a critical acceleration. As we explained in [38] to which the reader 
is referred for further details these limitations are due to the calculational scheme used 
here ("cutoff identification", etc.). Usually this scheme can provide a first qualitative or 
semi-quantitative understanding, but if one wants to go beyond this first approximation, 
a full fledged calculation of r[(7^,^] would be necessary which is well beyond our present 
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technical possibilities. 



6 Summary and conclusion 

In this paper we assumed that Quantum Einstein Gravity correctly describes gravity 
on all length scales and tried to identify the underlying RG trajectory realized in Nature. 
Along this trajectory, we found a regime where the running of the gravitational parameters 
is unmeasurably small, the domain of classical General Relativity. The renormalization 
effects become strong both at momentum scales k much larger (UV regime) and much 
smaller (IR regime) than those of the GR regime. In the UV regime, for k oo, the tra- 
jectory approaches a non-Gaussian fixed point. By definition, QEG is the "asymptotically 
safe" theory whose bare action equals the fixed point action, guaranteeing its nonpcrtur- 
bative renormalizability. We analyzed a potential IR divergence (instability) in the fiow 
equation which occurs in the Einstein-Hilbert truncation already, and we argued that it 
triggers significant renormalization effects in the IR, i. e. at distances larger than those 
where classical GR was successfully tested. In particular, Newton's constant increases for 
decreasing k so that G should appear to be larger at galactic or cosmological scales than 
in a terrestrial or solar system "laboratory" . Looking for possible manifestations of this 
effect we analyzed the possibility that the almost flat plateaus observed in the rotation 
curves of spiral galaxies, usually attributed to the presence of dark matter, are actually 
due to the RG running of Newton's constant. While the galactic regime is inaccessible 
with the Einstein-Hilbert truncation, a phenomenological analysis in the framework of the 
improved action approach revealed that a power law-type scale dependence G{k) oc k~'^, 
corresponding to an (approximately) constant anomalous dimension 77N, leads to rotation 
curves which are in qualitative agreement with those found in Nature. In particular they 
become perfectly flat at sufficiently large distances. 

As a by-product we arrived at a new understanding of the cosmological constant 
problem: By analyzing the RG trajectories among which Nature could have chosen we 
saw that there exist no trajectories which, at the same time, would predict a long classical 
regime and a large cosmological constant. If there is a Universe inhabited by human beings 
who manage to verify that gravity behaves according to classical General Relativity over 
many different length scales, they unavoidably will observe a cosmological constant which 
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is extremely small compared to the value of 1/G = mpj they measure. Conversely, if A 
is not small, it is also not constant, so GR is not valid, and there is no point wondering 
about the value of a non-constant "constant" . This argument does not seem to require a 
detailed understanding of the IR effects since the data indicate that in our world, between 
the end of the GR regime and the Hubble scale, G and A changed by comparatively small 
factors only. 

As for the IR effects, it will be exciting to see whether the properties of the RG 
trajectory of QEG which we predicted can be confirmed by more advanced ab initio 
calculations. Clearly more general truncations for the effective average action, involving 
nonlocal invariants for instance, should be explored for this purpose. In view of recent 
progress [18] made on causal (Lorentzian) dynamical triangulations [17] one also may hope 
that at some point it will be possible to make contact with numerical simulations. This 
could lead to an independent confirmation of the picture we have drawn in the present 
paper. 
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